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Abstract 

We study the solutions of the Yang-Baxter equation associated to nineteen vertex models 
invariant by the parity-time symmetry from the perspective of algebraic geometry. We 
determine the form of the algebraic curves constraining the respective Boltzmann weights 
and found that they possess a universal structure. This allows us to classify the integrable 
manifolds in four different families reproducing three known models besides uncovering a 
novel nineteen vertex model in a unified way. The introduction of the spectral parameter 
on the weights is made via the parameterization of the fundamental algebraic curve which 
is a conic. The diagonalization of the transfer matrix of the new vertex model and its 
thermodynamic limit properties are discussed. We point out a connection between the form 
of the main curve and the nature of the excitations of the corresponding spin-1 chains. 
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1 Introduction 



An important concept in the theory of soluble two-dimensional lattice models of statistical me- 
chanics is to embed the respective transfer matrices into a family of pairwise commuting operators 
P]. Let Tl ( oj[ , . . . , Un ) denote the transfer matrix on a chain of size L with Boltzmann weights 



u\, . . . , uJn . This approach requires that the transfer matrix fulfill the following property, 



0. 



(1) 



for arbitrary L and weights ujf , . . . , Un for i = 1,2. 

At first sight it appears one needs to verify an infinite number of constraints for the weights 
since Eq.([T]) should be valid for any values of L. This is however not the case since Baxter argued 
that a finite number of local conditions on the weights are sufficient to assure the commutativity 
among distinct transfer matrices for arbitrary L [2] . These conditions can be written as a matrix 
equation whose structure depend much on the family of lattice model under consideration. In what 
follows we shall discuss them for a relevant class of lattice systems denominated vertex models. 

Lets us for example consider a vertex model on a square lattice of size L x L. The respective 
statistical configurations sit on the edges of the lattice. In the simplest case the number of states 
living on the horizontal and vertical edges are the same and take values on a set of integers num- 
bers 1, 2, . . . , N. Furthermore, to each vertex configuration a, b, c, d — 1, . . . , N is assigned the 
Boltzmann weight C (a, c\b, d) as defined in Figure [U The sufficient condition for the commuta- 
tivity of two distinct transfer matrices associated to a given vertex model for arbitrary L is the 
celebrated Yang-Baxter equation [2]. Considering the notation of Figured! these set of functional 
relations can be written as follows, 



JV 

E 

7li72,73=l 
N 



R (ai, 7i|a 2 , 72) (71, ^1^3, 73) £ (2) (72, ^2(73, ^3) 
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71>72,73=1 



where R(a, c\b, d) are the elements of an invertible iV 2 x iV 2 matrix often called R— matrix. 
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Figure 1: Elementary configuration in a vertex model on a square lattice. 

In order to classify solvable vertex models with given statistical configurations one has to find 
the possible solutions of the corresponding Yang-Baxter equation. In principle, we can consider 
this problem by following the method devised by Baxter for two states vertex models [2]. We 
start by eliminating the matrix elements R(a,b\c, d) with the help of suitable subset of relations 
derived from Eq.fl2J). The remaining functional equations will then depend only on the weights 
C^'(a, b\c, d) and C^^a, b\c, d) which need to be decoupled by means of the technique of separation 
of variables. This step is essential to reveal us the algebraic invariants constraining the Boltzmann 
weights parameter space. The parameterization of such manifolds provides us the dependence of 
Boltzmann weights on spectral parameters which is for instance useful to formulate the algebraic 
Bethe ansatz 011112]. 

The implementation of all the steps described above, even for a given fixed vertex state con- 
figuration, is in general a tantalizing problem in mathematical physics. The concrete results are 
mostly concentrated on vertex models having two states per edge, see for instance [6j [7J [8]. In 
this case the non-trivial models that have been uncovered are mainly associated to the algebraic 
manifolds of the asymmetric six vertex model [9], the symmetric eight vertex model [TO] and the 
so-called free-fermions systems [UJ. The difficulties of the problem increase with the number 
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N of states due to the proliferation of the possible allowed weights that are ultimately fixed by 
different classes of functional relations. The merit of this approach is that it makes possible an 
unambiguous classification of integrable vertex models with a given statistical configuration from 
first principles. 

The purpose of this paper is to begin a study of the Baxter method for three state models 
satisfying the so-called ice-rule. The respective statistical configuration leads us to the total 
number of nineteen Boltzmann weights. Here we shall consider a relevant subclass of such models 
whose weights are invariant by the joint action of the parity and time (PT) reversal symmetry. For 
this family of nineteen vertex models we show that the dependence of the underlying algebraic 
curves on the respective weights are rather universal. The possible solvable vertex models are 
classified by the distinct branches of certain invariants values entering the definition of these 
manifolds. This general analysis provides us the classification of nineteen vertex models in four 
different families. The first three of them have already been discovered in the context of integrable 
spin-1 XXZ chain [12], the quantum inverse scattering of the Mikhailov-Shabat model [13] and the 
quantum algebra U q [SU{2)] at roots of unity [HJ US EE]. Interesting enough, our results reveal 
that these models have the same underlying algebraic background despite their rather distinct 
quantum group origin. To the best of our knowledge the fourth uncovered vertex model is new in 
the literature. 

We have organized this paper as follows. For sake of completeness we review in next section 
the main characteristics of the nineteen vertex models. In section [3] we consider the analysis of 
the functional relations coming from the Yang-Baxter equation. We develop a systematic way 
to solve such Yang-Baxter relations leading us to determine the algebraic curves fulfilled by the 
Boltzmann weights. It turns out that the principal algebraic curve is a conic involving three basic 
independent weights. The remaining amplitudes of the parameter space are remarkably resolved in 
terms of ratios of polynomials depending on such basic weights. This makes it possible to classify 
the VT invariant nineteen vertex models from a unified perspective. In section H] we discuss the 
parameterization of the Boltzmann weights in terms of a spectral parameter and the associated 
spin-1 Hamiltonians. We note that the geometric form of the fundamental algebraic curve is 
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directly related to the nature of the excitations of the corresponding spin chains. In section [5] we 
present the eigenvalues of the transfer matrix associated to the new nineteen vertex model and 
its respective Bethe ansatz equation. We investigate the bulk properties of this model providing 
additional support to the mentioned relationship among the curve geometry and the behaviour 
of the spin-1 chain excitations. Our conclusions are presented in section [6j In the appendices we 
summarize some technical details useful for the understanding of the main text. 



2 The nineteen vertex model 

In this section we shall review the main features of the nineteen vertex model. This lattice model 
has three states per edge which here will be denoted by a, b,c,d = 0, ±. The allowed statistical 
configurations compatible with the ice-rule a + b = c + d lead us to nineteen different Boltzmann 
weights. These weights are represented in Figure |2] where the respective subscripts emphasizes 
the non-null charge sectors a + b. 
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Figure 2: The vertex configurations of a general nineteen vertex model on a square lattice. 



4 



In this paper we will investigate the integrable manifolds of nineteen vertex models whose 
weights are VT invariant. This invariance relates some of the weights leading us to subclass 
of three-state vertex models with fourteen distinct Boltzmann weights. More precisely, the VT 
symmetry imposes the following relationship among the weights of Figure [21 

b± = 6±, / = /, e = d, e = d. (3) 

We begin by introducing a common notation for the Boltzmann weights of such family of 
nineteen vertex models. Taking into account the parameter subspace (J3J) one can represent the 
content of the jC^-operators by the following matrix, 

£(*) = a ?ei,i ® e lfl + 6? [ei.i <8> e 2 , 2 + e 2j2 <g> e ltl ] + f {i) [ei,i g> e 3 , 3 + e 3 , 3 ® e lfl ] 

+ b® [e 2 ,2 ® e 3i3 + e 3i3 ® e 2j2 ] + t/ w e 2)2 ® e 2j2 + al l) e 3)3 ® e 3j3 + /i (i) ei j3 ® e 3) i 

+ ^ W e 3 ,i <8> e 1)3 + c^ei )2 <g) e 2j i + c^e 2 ,i <8> ei )2 + c^e 2)3 ® e 3)2 + cL J) e 3>2 <g> e 2j3 
+ [ei, 2 ® e 3i2 + e 2)3 <g> e 2 ,i] + [e 2 ,i <g> e 2 , 3 + e 3 , 2 ® ei, 2 ] , for i = 1, 2, (4) 

where e Qi b denote 3x3 Weyl matrices. 

By substituting the £W— operators (J4j) in the Yang-Baxter equation (j2J) we are able to deter- 
mine the null matrix elements of the R— matrix. Under the mild assumption that the weight a± 
are in general distinct from b± and /® it is not difficult to see that the R— matrix has the same 
form of the operators. For latter convenience we shall therefore express the R— matrix as, 

R = a+ 0) ei,i <g> e x ,i + ft™ [e 1A <g> e 2 , 2 + e 2j2 ® e hl ] + f {0) [e M ® e 3>3 + e 3j3 <g> e lfl ] 

+ &_ [e 2)2 <g> e 3>3 + e 3>3 g> e 2 , 2 ] + # (0) e 2>2 <g> e 2j2 + aL 0) e 3i3 g> e 3>3 + ft, (0) e li3 g> e 3j i 

+ /i (0) e 3 ,i <g> ei j3 + c^ei,2 ® e 2j i + c+ ) e 2> i ® ei )2 + c_ e 2 , 3 (g> e 3i2 + c^e 3j2 ® e 2>3 

+ d {0) [ei, 2 <8> e 3i2 + e 2i3 <g> e 2 ,i] + d (0) [e 2 ,i ® e 2>3 + e 3 , 2 g> ei, 2 ] • (5) 

We emphasize here that we are interested to classify genuine nineteen vertex models and 
therefore all the weights a± , b± , c±, c± , d^, d^, and foW f or i = 0, 1, 2 are assumed 

to be non-null. We recall that a classification of solvable three-state vertex models has been 



5 



attempted before in the literature [IT]. In the work [IT], more stringent symmetry conditions for 
the weights were assumed besides the fact that some of them could be null. As a result, the only 
strict nineteen vertex model found was the standard Fateev-Zamolodchikov spin one model [T2] . 

3 The Yang-Baxter Equation 

The purpose of this section is to investigate the functional relations for the Boltzmann weights 
which are derived by substituting the £W— operator (T4]) and the i?-matrix (J5J) structures into the 
Yang-Baxter equation (j2J). We find that the functional relations can be classified in terms of the 
number of distinct triple product of weights. It turns out the minimum number of triple products 
is two while the maximum is five. In Table [T] we summarize the number of different equations 
having two, three, four and five types of triple product of weights. Clearly, both the number and 
the structure of the functional relations to be analyzed are far more involving than that associated 
to the N = 2 state vertex model satisfying the ice-rule Q. We shall therefore start our discussion 
by first solving the simplest relations containing two triple products. 



Number of triple products 


Number of equations 


two 


6 


three 


36 


four 


57 


five 


24 



Table 1: The dependence of the number of distinct functional relations on the respective number 
of triple products. 



^n this case we have an asymmetric six- vertex model and we end up with six non-trivial relations containing 
only three triple products. Their solution is easily found by imposing that two distinct 3x3 determinants made 
out of the six relations are null [B] . 
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3.1 Two terms relations 

The six equations possessing only two different types of triple product of weights are given by, 

cM 2) - cfi 1 ^ = 0, (6) 

d (0) J(1) C P) _ J(0) d (l)g(2) = Q) (7) 

cfjd)^) _ 4°) d (D J(2) = o. (8) 

We first note that the apparent difference between c± and c± can be gauge away by a transfor- 
mation preserving the Yang-Baxter equation. Without losing generality we can set, 

cf = cf, for % = 0,1,2. (9) 

By substituting the result ((9]) in Eqs.( l7]l8|) we conclude that the weight becomes proportional 

to d®, 

^ = % for 2 = 0,1,2, (10) 

where \l/ is our first invariant value. 

We shall now turn our attention to the relations involving three types of triple products. 

3.2 Three terms relations 

By using the solution (|9]fT0l) we find that the thirty six relations possessing three triple products 
are pairwise equivalent leading us to eighteen distinct functional equations. They can be classified 
in terms of three different groups of six equations given by, 



j;roup G± 





+ cfa { l ) bl ) 


- cfbfa^ 


= 0, 


bfc^bf 


+ ci 0) ai 1) ci 2) 


- ag><g>ag> 


= 0, 




+ 4 0) «i 1) i 2) 


- agbgg 


= o, 



(11) 

(12) 
(13) 
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jroup (j± 



jroup G± 



d (0) 6 (l) & g) + / (0) d (l) c (2) _ rf(0 ) /( l) a (2) = Q) 

/<P)d«&£> + d^bfcf - b^d^a® = 0, 

w^d^ + /W6i 1} 4 2) - ^4 0) / (1) 4 2) = o, 

&£>&£>dW + c m d W /( 2 ) - af /Wd« = 0, 
4°)6£)d< 2 ) + tiPdP / (2) - ofd* 1 ^ = 0, 

m^d^d® + 4°)6« - c f /wjg) = o, 



(14) 
(15) 
(16) 

(17) 
(18) 
(19) 



where the subscript ± in Eqs.( jTTllT9l) means that each of them splits into two independent func- 
tional relations. 

Altogether we have eighteen linear homogeneous equations but only eight Boltzmann weights 
a± \b± ,c± , d^ and are at our disposal to be eliminated. Therefore, we have a high degree of 
over-determination to overcome making the solution of Eqs.( fTTiri9l) far from being trivial. We shall 
start our analysis by considering the group G±* of equations. They are similar to the functional 
equations underlying the symmetric six-vertex model and can be easily handled. We first eliminate 
the weights af± , b± in terms of c± with the help of Eqs. (1111)121) . As a result we obtain, 



JO) 



(i) 

± 



'i 2) 



r (2) 



1 2 



bj»_ 
JO) 



J1)J2) 
C± C± 



(20) 
(21) 



By substituting the above results (I20ll2ip in Eq. (fl3l) we find two separable curves that involve 
the weights a^,6± and c± for i = 1,2. Their solution leads us to invariants typical of six- vertex 
models, 





2 




2 




w 

a± 


+ 


~bf 




~M 



a^bT 



A+, for 



1,2, 



(22) 
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where A± are free invariant parameters. 

(2) 

We now turn our attention to the relations associated to the group G± . We start by first 
eliminating the weights and from Eqs. ( 114)15]) . Because they can be isolated from equa- 
tions possessing different charge sectors we need to keep track of their explicit expressions. The 
compatibility of such distinct solutions will be implemented subsequently. The expressions for S°> 
and /(°) are, 



&£>a£ 2) - a£>b. 



(1)1,(2) 



JO) 



Jl) 



f(0) 



fWgf _ fcWfcg) 



J2) 



1 2 



(23) 



(24) 



(o) 



where as before they are given in terms of the common weights c± 

We then substitute Eqs. ( 123)24"]) in Eq. (ll6j) and by taking into account that the Boltzmann 
weights are non-null we find the following relation, 



+ [/ (1) ] 2 W 2 >4 2) - 



.(2) 
''± 



+ 



.(2) 



-*[dW] 2 ai 2) ^ 2) 4 2) =0. 

(25) 

By using the previous invariant ([22]) on the second term of Eq. (]25]) one is able to split the weights 
with index i — 1 from those labeled by i = 2. The solution of Eq.(]25]) leads us to the new 
constraints, 



•,(i) 



+ [/«] 2 -A ± 6i 1) / (1) 



J2) (2) 

bi'c T 



(26) 



where A± are once again constant parameters. However, due to the consistency of the right hand 
side of Eq. (126]) they are related by, 

A + A_ = # 2 . (27) 

(2) 

In order to complete the analysis of the group G± we still need to impose the compatibility 
between the two possibilities for S ^ and derived from relations (123)241) . The consistency for 



the weight is easily resolved by fixing the relation between the amplitudes c + and c . By 
way of contrast the compatibility for weight requires us to identify the expressions for the 
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charge indices ± of Eq.( )24|) . The result of such identification is, 



(2? 



Fortunately, Eq.( )28l) becomes separable once we take into account the right hand side of Eq. (I26H . 
This allows us to relate the products by^c^ with b^cy and as a result Eq.( l2"gj) can be solved by 
the method of separation of variables. The solution is, 



A+aW _ 7,(1) 



b_ — by a" 'ey — aye 



(2)„(2) _ „(2)„(2) 
+ 



/CD 



n, 



(29) 



where Q is a new invariant value. 

(3) 

Before proceeding with the analysis of group G± we should pause to discuss the results 
obtained so far. The main feature of the last invariant (1291) is that it relates weights with different 
charge index. We expect therefore that the invariants values A±, A + ,\l/ and Q should not be 
independent of each other. One way to unveil such constraints is to proceed as follows. By using 
the left hand side of Eqs. (l26ll29p we can eliminate the weights S 1 ^ and b_\ namely 

2 _ ..._o „s „s 



[d«] 



2 _ J_ 



+ [/«] 2 -A + 6jV (1) , 



A, 



Qf 



(D 



(30) 



We now substitute the results (150]) on the left hand side of Eq. ([2"oT) with the charge index — . After 
using the relation (j2"7|) we obtain the expression, 



+ A + ^ - A_A + ) b ( +> + (A^ - A_A+fitf - ^ 2 + Q 2 ^ 2 ) = 0. 



(31) 



We next consider analogous approach to the weights having index % = 2. We first eliminate 

(2) (2) 

the weights by and c_ with the help of the right hand side of Eqs.( 

A + &i 2) cL 2) 



-,(2) 



.,(2) 



(2) (2) 

ay'cy 



* ci 2) 



(32) 



~ a?> + nb^ ' 

By substituting Eq.( l32l) on the invariant connecting the weights a_, b^ and c_ given by Eq.(l22 
one easily obtains, 



<(2) 



A_A,^&: 



(2) 



,(2) 



Qb 



(2) 



,(2) 



2Qa ( r'b 



(2) .(2) 



Q 2 



.(2) 



0. (33) 
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The above relation can be further simplified by using in the last term of Eq. (1331) the expression 
for the invariant A + , see Eq. (!22l) . By performing such simplification we find, 

^ + A + ^ - A_A+) af + (A 2 - A_A+ft# - ^ 2 + Q 2 ty 2 ) bf = 0. (34) 

We now reached a point in which two different possibilities emerged. The first one consists in 

the assumption that the pair of weights /^} and {a+\ b^}} are considered linearly dependent. 

b m a (2) 

From Eqs.( l3T|3"4")) we see that this latter hypothesis implies in the identity jpy = which together 
with Eq. f)24p lead us to conclude that the weight vanishes. Here we stress that we are looking 
for genuine nineteen vertex models and therefore such solution is disregarded. Thus, we are left 
with the second possibility which is simply to set the coefficients of Eqs. (l31ll34l) to zero. This 
condition implies that the invariants A±, A + , \1/ and Q are constrained by the following equations, 

+ A+^ - A_A+ = 0, (35) 

A 2 + - A_A+fM - ^ 2 + tt 2 ^ 2 = 0. (36) 

(3) 

Let us now discuss the solution of the functional equations associated to the group G±. For 
this group we see that the weights at our disposal to be eliminated are a± , b± and c± . However, 



they have already been computed by means of Eqs.( flTllT3l and therefore our task consists to 
make the equations of group G±* compatible with those of group g£ . The first step to solve this 

(3) 

problem is to assure from the start that all the six functional relations of group G± are indeed 
satisfied. This is done by eliminating the weights af± and b±* from Eqs. (117fl8p . 



1 



dWd^ 1 6 W 6 P) _ /(l) f(2) 

b (0) 1 f &£)/(!) [rf(2)] 2 _ [rf(D] 2 4 2 )/(2) 



(37) 



(38) 



4°) d^d^ \ fo ( = 1 ) b ( t 2 )_ / (D / (2) 

By substituting the expressions f l3"Tf3"8"|) in the last equation of the group G± \ i.e. Eq.f lT9~j) . we 
obtain for the charge + sector, 



b (i) uw h m _ ^ ^ i _ h m /(2) | 6 w 6 « _ * [ rf d)] 2 | _ [ f w] 2 h {2) + b m /« [f^y = o, 

(39) 



n 



while for the charge — sector one finds, 

| 6 (2)^2) _ y J rf(2 )j2| _ ^ /(2) _ ^ ^ {1) j2 j _ ^(1)^(2) /(2) + & (1) j /(a) j2 = ^ 

(40) 

While Eqs. (I39|HU1) are not individually separable it turns out that suitable linear combinations of 
such equations can be split on the indices i = 1, 2. In fact, by adding and subtracting Eqs. (139(40]) 
we found that they are solved provided the following constraint is verified, 



9±, for i = l,2, (41) 



b w b w _ ^ [ d (0]- ± [f(i)y 

where Q± are additional invariants. 

Once again the invariant fj4"Tj) connects weights carrying distinct charge sectors ±. We have 
therefore to repeat the same analysis we did for the previous invariant Q, see Eq. (l29|) . As before, 
the expression fT4T]) for index % = 1 is simplified with the help of the weights and 6^ given by 
Eq. (l30p . This leads us to the following polynomial identity, 

[A+ ± * - ^0±(A + + Q)] + [±n% =f 0±(A+ t *)] = 0. (42) 

Considering that we are looking for linearly independent solutions for the weights b+ and 
we are required to set the coefficients of Eq.( ]42l to zero. In other words we have the additional 
restrictions, 

A+±^-^0±(A + + ft) = 0, (43) 

and 

±Of T 9±(A + Tf)=0. (44) 



It turns out that by solving Eq. ()43l for the invariants Q± and by substituting the result in 
the companion relation (j44l one finds that it is trivially satisfied once we consider the previous 
constraint fl35|) for the invariant A_. This means that Eqs. (j43|44p are both solved provided that 
we choose, 
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i2) (2) 

For the index % = 2 we once again have to use the expressions for the weights b_ and c_ 
given by Eq.( l32l . By substituting these weights in Eq. (HTT) we obtain two linear equations for the 

(2) 

weight a_ associated to the charge sectors ±. They are given by, 







\ 


[*e ± [d^] 2 + bff^ T e ± [f^] 2 } 









(46) 



The compatibility of such relations for the sectors ± and the fact that the weights should be 
not null lead us to the following relation among the weig hts bf\ and f( 2 \ 

2 



(e_ - e. 



l(2) 



+ [ 2 + M> (e+ + e_) [rf (2) ] +2(1- + 0_) fef/ (2) = o. 



(47) 



Taking into account the constraints ( )35|36|45p . one is able to show that Eq.( l47|) becomes equivalent 
to the following expression, 

2 



^(A+ + n) 



■,(2) 



The unique solution of Eq.()48l) that does not drive us to null weights 



(4£ 



',(2) 



+ [f^Y-A + b?fW 



[d( 2 



A. 



(49) 



Interesting enough, Eq. (|49|) generalizes the left hand side of the invariant (|26|) to the weights 
with index i = 2. At this point we observe this fact also works for the invariant A_ involving the 
weights b^ 2 \ d^ and f^ 2 \ In fact, by using the expression for the weights b^ 2 \ cl\ ol^ coming 



from Eqs. (132|46j) besides the constraints f l2Tj35f36|l45]) among the invariants as well as Eq.fl49J) to 
eliminate the weight [d^ 2 ^] 2 one is able to verify that, 



<(2) 



+ [/( 2 )] 2 -A_&L 2) / (2) 



[d^x 



A. 



(50) 



(3) . 

We can now return to discuss the solution of the remaining equations of group G±. At this 
stage we have just to match the weights a!± and b± obtained from distinct pair of relations 



2 Note that * = implies = while A+ + il = is equivalent to ^ 2 V 2) - * [d {2) ] 2 ± [/ (2) ] 2 = 0. The 
latter identity implies that = 0. 
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associated to the groups G± and G*± . In other words, we have to make Eq.( l20l) compatible 
with Eq.(l371) as well as Eq.(l2Tl) consistent with Eq. (l38l) . This require us to solve the following 
functional relations, 



b^af-a^' 



b £) b m_ f m f ® 



J1)J2) 



6«/ (1) ^ (2) ] 2 -^ (1) ] 2 4 2) / (2) }, (51) 



l(2) 



J2) 



(1) f (2) 



f W f 



[rf (2) ] 2 -[rf (1) ] 2 [/ {2) ] ; 



(52) 



Further progress are made by squaring both side of the expressions (j51|52p . This operation makes 
it possible to use the invariants A± (1221) to eliminate the weights c± as well as Eqs. (I26|l4ljf5"0]) 
to extract [<i^] 2 for both i — 1,2. By performing such two step procedure we are able to show 
that Eqs. (151l52j) become proportional to the expression, 



—A. 



ai 2) + 



4 2) / (2) 



2 



(!) f (1) 



- a£>f 



(!) f (1) 



a£ J f 



l(2) 



0. 



l(2) 



-«i 2) / (2) 



+ 



(53) 



By inspecting Eq.fl53l we conclude that such relation can indeed be separated leading us to our 
last invariant associated to the three terms relations, 

a&h£ ) -A ± agf® + b , £f® 



bf 



r± for i = 1,2. 



(54) 



where r± are additional variables. 

The final step of our analysis consists in matching the different charge sectors ± of the invariant 
(1541) . This task involves the manipulation of cumbersome expressions and the respective technical 
details are presented in Appendix A. The condition of consistency turns out to be a constraint 
among the invariants r±, A + and Q whose expression is rather simple, namely 



r_A+ - #(r+ + n) = o. 



(55) 
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In addition, an important byproduct of the analysis performed in Appendix A is that the only 
independent weights are a$ , b^} and c+ . In fact, remaining amplitudes entering the three terms 
relations can be written in terms of such weights. As a result, the weights a®, b_ , c^}, g?(0 and 
/(0 are given in terms of ratios of polynomials whose degrees are at most two. In what follows we 
shall present such relations since they are going to be useful later on. Following Appendix A, the 
simplified expressions for weights al , b_, c_ are, 

^ 2 (a? + n&J ) )[(A + - r + + n)a { ; ] - (i + r+tyb*?] 



.0 #.(0i 



',(0 



A^[(A + -r + )a^-^ j 

*[(A + - r + + 0)0? - (l + r + o)#] 



A+[(A + -rjo2 ; -6: 



,(*) 

(0 M i r rmCOl 



* 2 [(A + - r + + 0)0? - (l + r + o)6^] (i) 



A 2 + [(A- 



T i )a 



(0 



while for weights and /(0 we have 



i-A,r,+r 2 



M„(0 



A, 



A, + T + ) a? + 



/ 



(0 



+ ' 



(56) 
(57) 
(58) 



« 



l(0 



(A- 



l(0 ' 



(59) 



We would like to conclude this section with the following comments. First it is important to 
stress that all the invariants obtained so far are also valid for index i = 0. This result is verified 
by using the expressions of the eliminated weights a± , b± , c± , d^ and /(°) on the form of the 
respective invariants. As a consequence of that, the above expressions (I56H59I) remain valid for 
the weights cS®\ cl\ d^ and f(°\ We next note that out of ten possible invariants values 
A±, A±, T±, 0±, Q and \I> we end up with only four free parameters because of the six constraints 
(I27)35)36|45)55p . For sake of completeness we summarized such conclusions on Figure [3j 



3.3 Four terms relations 



After using the two terms solution 



the total number of four terms relations is reduced to 



fifty one equations. The majority of these relations depend on the weights and hS 1 ' that 

are still to be determined. The only exception is a factorizable functional equation given by, 



o. 



(60) 
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,«] +[/W] 2 - A ± ^/W ft W c W 



n 



6 CO ±6 (0] f(0 



e 



± 



Figure 3: The invariant structures solving the three terms functional relations f lTTTfT9|) valid for 
i = 0, 1, 2 whose values are constrained by Eqs. (I27|35ll36ll45|55p . 

The other fifty equations can be classified in five different groups characterized by the number 
and the type of unknowns weights hS l \ present in such functional relations. We have 
equations involving only one of the weights, the conjugate pair of amplitudes and as well 
as all the weights g( % \ together. In what follows we list these distinct groups of relations, 



group G± 



group G± 



- c^c^b^ + HfdF>S»b%> + g<%g<§> - bfgVcf = 0, (61) 

- b^c^c^ + M&V 1 ^ - cfgUb® + (4 0) 6iV J = 0, (62) 

- A«c£ 2) + cfb^d^ - g^S 1 ^ + b^d^gW = 0, (63) 

h (o) {b w b m _ 6 a) 6 W) + d (o) d v {( fi) _ C W) = o, (64) 
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h^b 1 ® - b {0 ¥ 2) ) - c^cVc® + cfc^c? = 0, (65) 



group G±* 



(7) 



• group G± 



fcW( 6 W 6 W _ 6 <_°> 6 «) + d^S 2 \cf - c^) = 0, (66) 

- d^S%^> + <£<V>b%> - h^cf + bfhF><%> = 0, (67) 

cfh^bf + ft^c^ - bfd^d^ - cfb^hM = 0, (68) 

a^h^a^ - dWgdP - f {0) h^ - hWaghP = 0, (69) 

h^ib^ - b<l¥?) + tfW^ci? - c?>) = 0, (70) 

fc(D( 6 W 6 W - 6 W 6 <2>) + c (0) c (l) c (2) _ C W C (1) C (2) = Q; 

fcW(_ 6 W 6 W + 6 <_°> 6 «) + ^V 1 )^^ - cf) = 0, (72) 

tfW% 2 > - cfc^ - blWtg + = 0, (73) 

_ 4°>/L(D 6 W - ^c^c^ + ^dW^ + c^/^ = 0, (74) 

- afh^af + /W/iW /( 2 ) + &dP><gdP> + /^a^ = 0, (75) 

_ 1)d (0) c W d (2) + ^(0)^(1)^(2) _ ^(0)^(1)^(2) = Qj (76) 

^(o)dd) 6 w + d(o) 6 «4 2 ) _ *<4°>&«d( 2 > - fcWjd)^) = 0> (77) 

_ fc(o)d(i) 6 W _ M/d(o)6«ci 2) + c^'d* 2 ) + mfdF>hV> = 0, (78) 

fr<0) /(D a g) _ h (o) a W ;(2) _ ^rf(o) c « d (2) _ / (0) /l (D/ i (2) = 0j (79) 

_ ~ h (0) f (l) a (2) + ~ h{ 0) a (l) /( 2) + ^(0)^(2) + / (0)/ i (l)^(2) = 0> (8Q) 

h (°)~ h m /(2) + *d(0) c « d (2) + / (0) a «/ l (2) _ a (0) /( l)£( 2 ) = Qj (gl) 

_ ftW^l) ; (2) _ ^rf(0) c « d (2) _ f (0) a W h m + fl (0) /( l) ft (2) = Qj (g2) 
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/-y(8) 

jroup (j± 



c w rf (D4 2 ) _ h (°) a u d w - g<n = o, 

a (f> d (D c g) _ 5 (0) c « d (2) _ tfdWfcW ;(2) _ rf (o) a a)^(2) = Qj 

«<!> /W _ Vc^dPtg + VgMfgdP) + VdPcghW = 0, 

h (o )c (i) /(2) _ c (o) /( i )c (2) + ^ d (o) p( i) d (2) + f (o) c m~ h m = 0j 

4 0) /«4 2) - tfdVW 2 ) - / (2) - /^i 1 ^ = o. 



J3) 

14) 
15) 
16) 
J7) 



From Eq. (160j) we see that we have to deal with at least two possible branches. Either \l/ 2 = 1 or 
c+ = c_ since the weights are assumed non-null. This information, however, is not necessary 
to obtain the invariant values associated to the weights gW. Indeed, we observe that Eqs. (16 If 621) 
from group G { ± ] contain the triple products c£ 1} c£ 1} , c^c^fe^, ^fd^d^b^ and tf&^rfMd* 2 ) 
which are also present in the previously solved three terms relations, see Eqs. (lllfl3fl6fl9p . By 
eliminating the first triple products of the latter mentioned relations and substituting them in 
Eqs. ()61f62l) we obtain the following expressions, 



■,(o) 



l± 



4 0) 



,(o) 



6 W _ f (o) b W + 6 (0) /( i) = 0j 



and 



4 1} - 9 {l) 



4 2) - g {2) 



+ f O)0 - fig) /<*) = o. 



i9) 



(90) 



An essential step to separate Eqs. (I89f90l) is to use the invariant f2 which links the weights b^} and 



6_ . Indeed, from the third box of Figure [3] one is able to eliminate b± in terms of b$ and 







r(i) 



By substituting the result in Eqs. (I89II90I) and after few algebraic manipulations which includes the 
use of identity (1270 we find that the invariant associated to the weight g® is, 



-9 {i) + off + 



for z = 0,1,2, 



(91) 



where are constant parameters. 



Our next task is to assure the consistency of Eq. (l9TT) as far as the charge sectors ± are con- 
cerned. This requires us to impose that the expressions for g( % > coming from the different charge 
sectors are the same, namely 



(i) 

ay — a 



»+(*__*_) /(0_Ato) 6 W +A to) 6 (0 



for i = 0,1,2. 



(92) 



In order to solve Eq.(l92l) we substitute the expressions for the weights ay, by and derived in 
previous section, see Eqs.f l56|l57f59|) . After some cumbersome simplifications we find that Eq.f l92|) 
becomes a polynomial relation on the amplitudes a+ and having the following formf , 



Ax 



+ A 2 a { l ) b { ! ) +A 



for i = 1,2. 



(93) 



The coefficients Ax, A 2 and A% depend solely on some combinations of certain invariants and 
should vanish to assure the validity of Eq. (|92|) . Their expressions after using the constraint (I2T1) 
are, 

a 1 = (A + -r + )A^-(A + -r + + ^)^ 2 = o, 

A 2 = -A 3 + - [1 + A? ( A + - r + )] A 2 + ^ + [1 + A (9) (A + - T + + ft)] A + ^ 2 
+ [l-fi(A + -2r+ + fi)]^ 3 = 0, 

a 3 = A^^ + r+A 3 -[A (3) + r + + A (9) r + fi]A + ^ 2 + (i + r + fi)fi^ 3 = o. (94) 

The next step consists in fixing the form of invariants associated to the Boltzmann weights and 
hS % \ To make progress in this direction we have to use explicitly the two branches data encoded 
in Eq.(|60|). 



3.3.1 Branch 1 

This branch is chosen by setting ^ = ±1. This constraint allows us to eliminate the weight 
from Eq. fl76|) 

h<X)h&) ... 

(95) 



h (0) = y_^ h (o)_ 



3 For i = the compatibility equation (|92|) is automatically satisfied once A\ = = A3 = and provided that 
we take into account the relations (|27|35|36[) among the invariants. 
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By substituting the result ( 195|) in the pair of equations (I69II75P and (I67f73p we find that their linear 
combination fix a relation between m 1 ' and h^' which is, 

fcW = h {l) for i = 0,1,2. (96) 

As a consequence of Eq. (1961) the equations of the group G± and become the same as well as 
the number of independent relations of groups G±> and G v ±> are reduced to half. At this point we 
can apply the same strategy used to solve the weight g®. We first eliminate the triple products 
of weights d^b±c± and c^b^d^ from Eqs. fl!5|18j) and by substituting the results in Eq. ( 177)) 
we find a separable expression for and , namely 

6f) " fcf) - A± ' (97) 

where are constants of separability. 

Once again we have to implement the compatibility among relations associated to different 
charge sectors ±. This matching is performed in the same way we did for the weight As a 
result we find the and derived from the charge sectors ± agree provided that the following 
relations between invariants are satisfied, 

= + r + Q)(A {h) A + y - Q), (98) 
A-+ 

fi = (A + -r + )(A 2 + -l), (99) 

(A 2 + - 1){1 + (A + - T + )[(A + - r + )(A^ - 1) - A W A + M/]} = 0. (100) 

From Eq. (11001) we observe that branch 1 splits in two distinct families since this relation admits 
two possible solutions. We shall denominate such possibilities branches 1A and IB. For branch 
1A the value of A + is fixed by, 

A+ = 1, (101) 
while branch IB is defined by solving Eq. fllOOp for 

A + *(A + -r + ) • {W2} 
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It turns out that for both branches we are able to manipulate Eqs. (167(68]) in order to determine 
the only remaining weight It is given by an expression similar to that found for and 
which is, 

-m — = aJ ) . (103) 

We now have reached a point in which all the Boltzmann weights have been determined in 
terms of the invariants values and the weights a+ , and c+ . By substituting them in the four 
terms relations not used so far and using the invariant (122]) to eliminate the weight [c^] 2 we obtain 
polynomial equations depending on variables a® and b+ whose coefficients are functions of the 
invariants values. By setting these coefficients to zero and by considering the previous constraints 
(I27ll35ll36|45|55p as well as Eqs. (l94|l9~8|99|ll00p we are able to compute the invariants values. It 
turns out that once Eqs. (I63f83l) are satisfied all the other remaining relations involving four terms 
are automatically fulfilled. The final result of such analysis is summarized in Table [2J 

Here we remark that the cases \1/ = — 1 or A + = — 1 leads us to specific invariants values which 
can be reproduced in the context of solution IB and another branch that are going to be discussed 
in next subsection. The technical details concerning such special situations have been collected in 
Appendix B. We finally observe that branches 1A and IB have a unique free parameter chosen to 
be A+. However, only the branch 1A is invariant under the charge symmetry +-<->■—. 

3.3.2 Branch 2 

Another possibility to satisfy Eq. fl60|) is to set This condition together with the previous 

relation (|58p lead us to a polynomial equation for the weights a+ and , 

[(A 2 + - ^ 2 )(A+ - r + ) - m 2 ]a { ^ + \-A 2 + + (1 + T+fi)^ 2 ]^ = O. (104) 

The linear combination (11041) is fulfilled for arbitrary and 6+ by imposing that its coefficients 
are null. This fixes the values of the invariants Q and A^ 

Q = and A + = (105) 



4 Note that the possible solution A + = T + + T, 1 is disregarded since it leads us to = 0, see Eg . (|55)) . 
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Invariants 


Branch 1A 


Branch IB 


A + 


free 


free 


A_ 


A+ 


-A++e 1 v / 3 v /4-A^ 
2 


A + 


1 


2e lv / 3 v /4-A2 
3A++e lv / 3 v /4-A2_ 


A_ 


1 


>/3A++ei v /4-A2 
2ei v /4-A2_ 




1 


1 







6-3A^-ei\/3A +v /4-A^ 
3A++eiv / 3 v /4-A2_ 


r + 


A+ + ei 


3A+-e lv / 3v / 4-A2_ 
6 


r_ 


A+ + ei 


-3A + +e lv / 3 v /4-A2_ 
6 


©+ 


2 
AT 


A + +eiv / 3 v /4-A2 
2 


e_ 





-3A++ei V / 3 A /4-A2_ 
6 


A? 


A+-ei 


V4-A2 


A (9) 


A+-ei 


(4-A2)(A ++£1 ^3 v /4-A^) 


eiv / 3A +v /4-A2+4-A2 




A + 


A + 




A + 


2ei 



Table 2: The invariant values for branch 1 where the discrete variable e\ = ±1. 
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By substituting fll(J5j) back in the expressions for the weights o_ , b_ and c_ given in Eqs.( l56ll571l58l) 
we find the rather simple relations, 

a W 6 W c W 

— = 1, - 7y = ±, — = 1 for z = 0,l,2. (106) 

aY o+ cY 

An immediate consequence of Eq. (jl06j) is that we can implement several simplifications on the 
four terms equations. First, we clearly see that Eqs. (l64(65|l66|l70|l71|l72p are automatically satisfied, 
besides that the number of independent equations reduces dramatically. For groups Gi 5) and Gi 6) 

(7) (8) 

we just have three distinct relations while for G± and G± the number of equations is reduced 
to half. We are now in position to determine the weights and hS l > by using the same method 
that fixed In fact, by eliminating the triple products c+ c+ , c+*c+b+\ tyd^d^b^ and 
m^d^d^ with the help of Eqs. (|ll|ll3fl6fl9ft and by substitut ing them in Eqs. (167(681) we are 
able to obtain separable functional relations for the weight /&W . The solution of such relations are, 

_ tfr^W + + f(i) ... 

+ +/ = Af J for i = 0, 1, 2. (107) 

By the same token the weight /i 1 -*- 1 can also be calculated. By using the above explained 
approach but now for Eqs. (173|74p one finds, 

-W0 +a » + »/» =A$) for ! = _ ii2 ^ 



At this point all the Boltzmann weights for branch 2 have been determined. The remaining 
task is to substitute them in Eqs. (I69|75|) as well in the relations of groups G± and G± to fix 
the invariants values. Once again we have to consider two possible branches since from the very 
beginning Eq. (j!05j) has two allowed solutions. Taking into account all previous constraints among 
the invariants our final results are summarized in Table [3j 



23 



branch 


Branch 2A 


Branch 2B 


A+ 


free 


eiV3 


A_ 


A + 


-ei\/3 


A+ 


2-A^+e 2 A +1i /A2 -4 


u; 


2 




2-A2_+<e 2 A +a /A2 -4 


u; 


2 




2-A^+e 2 A +A /A2 -4 


—a; 


2 


n 








r+ 


ei 


ei 


r_ 


ei 




©+ 


2 
A+ 





0_ 







A? 


A+-ei 







A+-61 





A? 








Af 









Table 3: The invariant values for branch 2 where u = exp (^^p) and the discrete variables e x = 
e 2 = ±. 



Note that only branch 2A has a free parameter and is invariant under charge conjugation. We 
would like to conclude this section with the following remark. Although the algebraic curves for 
the weie hts and h® come from different equations for branches 1 and 2, their final expressions 
can be put in a unified form. In fact, they can be written as, 

+ + /«■) _ a£)> -W + Jg + vf® = A < s , for . = ^ 2 



6» * ' 6« 



The general result (I109p recovers branch 1 by simply making the identification A^_ = A 



(h) _ A (h) 



± 



while for branch 2 we have to consider Eq. (11061) and therefore we have to impose = = ±. 
In Table H] we summarized the invariants obtained in the analysis of the four terms functional 
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relations. 





* J- 


















&« 








-fcW+a^ +*/(*) 
6<*> 


= A w 



Figure 4: The invariant structure associated with the four terms equations valid for i — 0, 1,2. 
The expressions for the invariants A^ , A±^ and A± ^ are given in Tables [2] and [31 

3.4 Five terms relations 

The number of the five terms functional equations remains unchanged after the solution of the 
two terms relations. Below we present the structure of the corresponding twenty four equations, 

4 0) 4M 2) - d®h<»dF> - 0<°><4V° + <4% (1) i 2) - &d<® = o (no) 

4°) C «4 2 ) - ^d^h^d® - gMcggW + c^W^ - ^a?^ = (HI) 

m^b^d® + c^d^hW - d^c^gW -Vf^a^dW - h^h^d^ = (112) 

b { M ] d (2) + *<g>dF>hF> - VdWfggW - f^a^d® - m^k^d® = (113) 

tfdWft^ftg) - H> A« - gWtgdP - ^fcW/iM + fcW^cgi = (114) 

dC0) 6 W 6 g) _ d (o) a a) /( 2) _ % (0) c W d (2) _ ^(o) ^(i)^(2) + ^^(o) d (D c g) = Q (n5) 

6 £°) c « fe £) + ^/iWcg) - /^c* 1 ' /( 2 ) - fcW^/jM - df®gV<F> = (116) 

b (V c M b W + c f~ h (i) c (2) _ /( o) c W /( 2) _ ^(0)^1)^(2) _ ^2^(0)^(1)^(2) = Q (n7) 
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X^(0) 6 « fe (2) + ^(0) d (l) c (2) _ f (0) c (l) d (2) _ m (0) g{ l) g (2) _ y~ h (°) c U d (2) = Q (llg) 
d (%U b W + g(°) d W^ - * jWcW^ 2 ) - W _ hF><V>dF> = (119) 

bfb^d^ - /M + cfd^gW - gWgMjn - e^cgV 2 ) = (120) 

mbfb { l ] d^ - + v^M^) _ H/^yV) - vI/^cW^) = o (121) 

The above relations do not involve any new Boltzmann weights. As a consequence of that, 
we can substitute the weights expressions obtained in the two previous sections in Eqs. (lllOH121l) 
and search for further restrictions on the invariants values for branches 1 and 2. Remarkably 
enough, after the three and four terms relations are solved all the above five terms equations are 
automatically satisfied. Therefore, the four families of vertex models defined by the invariants 
given in Tables [2] and [3] are exactly integrable. 



4 Parameterization and Hamiltonian 

We shall present here the parameterization of the weights associated to the four distinct integrable 
manifolds of previous section. As discussed in section [3] we have only one fundamental algebraic 
curve involving the weights a!+ , b^} and c+ . The corresponding curve (1221) can be rewritten in 
the form of a conic, 

* - ^rVi) -(^r-A Vi = !. ( 122 ) 



2 7 V 4 

where the new variables %i and yi are related to the weights by, 

Xi = "Tip y* = li) • ( 123 ) 

c\' c\' 

The rational parameterization of conies involves only one spectral parameter Aj and can be done 

through hyperbolic functions, namely 

A + sinh(Aj) 
x i ~ -TTVi = cosh(Ai), y t = - (124) 

2 

For later convenience we introduce the following definition for the invariant A+, 

A+ = 2cosh(7). (125) 
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By substituting Eq. (ll23j) in Eq. (I124j) and after using the definition (I125[) one finds, 

„«_i ,(») _ sinh(Aj) (i) _ sinh( 7 ) 

+ ~ ' + ~ sinh (A, + 7) ' + " sinh(A, + 7) ' l ^ bj 

where a!+ has been fixed by freedom of an overall normalization. 

Before proceeding we stress that under the above parameterization the Yang-Baxter equation 
is additive as far as the spectral parameters Aj are concerned, 

A = A x - A 2 . (127) 

The remaining weights can be determinate in terms of a+ , b+ and c+ with the help of Eqs. 
( I56|57|58|59|91|109p . In what follows we list the simplified expressions of the weights associated 
to the four different integrable vertex models defined by the invariants values given in Tables [2] 
and El 

• Branch 1A 

This vertex model falls in the family of the factorized spin-1 S— matrix found by Zamolodchikov 
and Fateev [12] with the additional presence of the discrete variables e\ = ±1. We remark that 
the transfer matrix spectrum for t\ = +1 can be related to that with e\ — —1 only for L even. 
The respective weights are given by, 



,(») _ sinh(Aj) ( j) _ sinh( 7 ) 

smh (Ai + 7) smh(Aj + 7) 



= d (i) = ± sinh (7) smh(Aj) (i) = sinh(A, -J + j) sinh(Aj) 

sinh(Aj + | + 7) sinh(Aj + 7) ' sinh(Aj + \ + 7) sinh(Aj + 7) ' 

w _ -2 coshQ - 7) + cosh(f + 7) + cosh(2A 8 + j - 7) , , 

5 2sinh(Ai + f + 7)sinh(Ai + 7 ) ' 1 ] 

W = fcW = 2cosh(j-7)sinh 2 (| + 7) (131) 
sinh(Ai + f + 7) sinh(Ai + 7) ' 

where 7 = ^(1 - ei). 
• Branch IB 
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The weights of this branch were obtained on the context of the representation of the U q [SU(2)] 

quantum algebra when q is at roots of unity [TH [T5| [T6] . Note that this model is not charge 

invariant and the expressions for the weights are, 

(0 , _ sinh(Aj) a) sinh (7) 

a+ ~ i; &+ " sinh (A; + 7) ' + " sinh(A, + 7) ' [ ' 

a (i) _ sinh(Aj -7 + 70) sinh(Aj. - 7) _ sinh(7 - Aj sinh(Aj) 

sinh(Aj + 7 - 70) sinh(Aj + 7) ' " sinh(Ai + 7 - 70) sinh(Ai + 7) ' 



c (<) = sinh(7o - 7) sinh(A t - 7) = ^ = ± y/smh(j) sinh( 7 - 7o ) sinh(A^) 

sinh(Aj + 7 - 70) sinh(Aj + 7) ' sinh(Aj + 7 - 70) sinh(Ai + 7) ' 

_ sinh(Aj - 70 ) sinh(Aj) (i) _ -1 + cosh(2A^ + 70) + cosh(27 - 70) 
sinh(Aj + 7-70) sinh(Aj + 7) ' 2 sinh(Aj + 7-70) sinh(Aj + 7) 

ftO = h ® = sinh(7)sinh(7 - 7 ) (136) 
srnh(A, + 7-70) smh(Ai + 7) 

where 70 = yei. 
• Branch 2A 

Apart from the presence of the extra discrete variables t\ = e 2 = ±1 this vertex model is equivalent 

to the so-called Izergin-Korepin R— matrix [13]. The transfer matrix spectrum of this model is 

however independent of the parameter e 2 . The explicit form of the weights are, 

(i) _ 1 Ai) _ sinh(Aj) w _ sinh (7) 

smh (A, + 7) srnh(Aj + 7 J 



d>l) = ^UirM^M & = _ exp ( , ,(,, 
cosh(Aj + -f + 7) smh(Ai + 7) 



:i38) 



cosh(Aj + % + 7) sinh(Aj 
cosh(Aj + ^ + 7) sinh(Aj + 7) 



/ W = , 3-v 2 rv^TTT 7~~\ i ( 139 ) 

^(i) = ~ Sinh(| + 7) - omnyy - ,j t q^lx^j- t /; -r -p j- - /; 



(i) _ cosh(Aj + ^ + 7) sinh(Aj. + 7) - exp (2e 2 7) cosh(A^ + | + 7) sinh(A^) 

~ cosh(A, + f + 7) sinh(Ai + 7) ' ' 



~ (i) _ cosh(Ai + 4? + 7) sinh(Aj + 7) - exp (-2e 27 ) cosh(Ai + \ +7) sinh(A;j 

~ cosh(A i + f + 7) sinh(Ai + 7) ' { '"~ 

where as before 7 = ^(1 — ei). 
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• Branch 2B 



The nineteen vertex model defined by the weights given below appears to be new in the literature. 
This model violates charge conjugation symmetry since = —b^\}. The respective weights are, 



#=1, *4°=± .ty v 4°= .ty v ( i4s ) 

± smh(Ai + 7) ± smh(Ai + 7) 



^^ r^rl t - ^» = — p(^)^'. ("4) 

2 cosh(A — 27) cosh(A) \ 3 / 

(i) = sinh(A i + 2 7 ) sinh(A^) (i) = cosh(2A t ) 45 

cosh(Ai — 27) cosh(Aj) ' 2 cosh(Aj — 27) cosh(Aj) ' 

h(i) = cosh(A, - 2 7 ) cosh(A i ) + exp (^p) sinh(A, + 2 7 ) sinh(A,) 

cosh(Aj — 27) cosh(Aj) 

_ cosh(A t - 2 7 ) cosh(A,) + exp (gg) sinh(A, + 2 7 ) sinh(A t ) 
cosh(Aj — 27) cosh(Aj) 

where 7 = y(l — ei) + y£i such that ei = e 2 = ±1. 

For sake of completeness we shall also present the expressions of the one-dimensional spin 
chains associated to these vertex models. At this point we observe that the matrix representation 
of the weights (jlj) at A = is exactly the operator of permutations on C 3 <E> C 3 . This property 
tells us that the respective Hamiltonians are given as a sum of next-neighbor terms, 

L 

H=J J2Hi, i+ i, (148) 

i=i 

where Jo is an arbitrary normalization. 

The two-body term H ifi+ i is obtained by the action of the permutator on the derivative of the 
respective D- 1 '— operators at \ = 0. Its general expression in terms of spin-1 matrices is, 

Hi ji+ i = J\SfS i+1 + JiSj Sf +1 + J 2 (SfS*) (S i+1 S? +1 ) + J2 (S* Si ) (Si +1 Sf +1 ) 
+ ^3 (St) {S i+1 ) + h (S { ) (S^i) + Ji (SfSlS i+l ) + J 4 (S-Si S^) 

+ j 5 (s?s- +1 s* +1 ) + J 5 (srsf +1 s+ +1 ) + h (st + Sf +1 ) + h 2 [(Sf) 2 + (S* +1 ) 2 } 

+ 5 1 S^St +1 + 6 2 (S^SU 1 ) 2 + h(Sn 2 SUi + ^SnSUi) 2 (149) 
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where S and S z are spin-1 matrices satisfying the SU(2) algebra. 

We note that in Eq. ( 1149 j) we have omitted any chemical potential terms that are canceled under 
periodic boundary conditions. In Table H] we show the values of the couplings entering Eq. ( 11491) 
for all the four branches. At this point we remark that branch 2B Hamiltonian is a particular case 
of an integrable spin-1 chain discovered in [18J. This equivalence occurs when the parameters 
and t p defined in [18] take the values t 3 = — 1 and t p = 

We conclude this section with the following comment. Some of the physical properties of the 
first three one-dimensional spin chains have already been investigated in the literature, see for 
instance [191 EU 1211 [22]. This knowlodge includes the nature of their excitation spectrum over 
the anti-ferromagnetic ground state. Considering these results we conclude that the spectrum is 
massless when the free parameter A + belongs to the interval — 2 < A + < 2. By way of contrast, 
outside this region the spectrum is expected to have a non zero gap whose value depends on A + . 
Remarkably enough, such change of physical behaviour is directly related to the variation of the 
geometrical form of the fundamental algebraic curve (11221) . In fact, in the massless regime the 
corresponding curve is closed and it is described by an ellipse while in the massive region the form 
of the respective curve is that of an open hyperbola. This relationship between the geometric 
form of the principal algebraic curve and the nature of the excitations suggests that the anti- 
ferromagnetic regime of the fourth spin-1 chain should be gapless since in this case A + = ±V3. 
In next section we shall present evidences supporting such prediction. 

5 Exact solution of Branch 2B 

We shall present here the exact diagonalization of the transfer matrix associated to weights of the 
novel vertex model defined as branch 2B. As usual the corresponding row-to-row transfer matrix 
T(A) can be written as the trace of an ordered product of £W — operators parameterized by the 
spectral parameter A. Recall that the expression for such operators in the terms of Weyl matrices 
is given by Eq.(T4]) while the corresponding weig hts a?,^,^, S % \ S l \ /W, and are 

given in Eqs. fll43H14T)) . 
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branch 


1A 


IB 


2A 


2B 


Ji 


± 1 


2s/r 


-4- 1 




2sinh(^+7) 


^~ 2cxp(e 2 7) cosh(^+7) 


~T 2 


L 
• J l 


L 

"1 


•Jl 




7, 
"1 


Jo 
u 2 


1 =F 1 


V / 3sinh(7-^-) 1 


1 _|_ sinh(e27) 


-r- ei£2 
"T 2 


sinh(7) ' sinh( ; 2-+7) 


2r "F y/f 


sinh(7) cos h(^+7) 


J 

J2 


7 

<h 


J 

<J2 


7 

^2 


7 

J2 


7 


ei 


i\/3e\ 


cosh(^+7) 


iV3ei 


4sinh(7) 


8r 


4 sinh(7) cosh( ^ +7) 


4 


T 

■h 


J3 




J3 


^3 


■h 


■h 


1 1 1 


1 _i_ exp(— €27) 




t 

2sinh(7-7o) 2yfr 


2smh( 7 ) ± 2cosh(^l+7) 


—i =F 


■h 


J'2 

2 




•h 


J4 


■h 


J2 
2 


1 1 

2sinh(7) ' 2y/r 




2 

-i ± 


■h 


J 2 
2 




-Ja 




Si 


2 cosh(7)+ei 





cosh(2e27) 


i\/3ei 


2sinh(7) 


2[-ei sinh(7)+sinh(27)] 


4 


S 2 


Si 





4ei cosh(7)— cosh(2t27) 
2[-ei sinh(7)+sinh(27)] 


3<Ji 


S3 








sinh(2e27) £ 
cosh(2e 2 7) 1 




Sa 








sinh(2e27) ^ 
cosh(2t27) ^ 


— ie 2 "\/35i 


hi 





sinh(27-7o) 
2r 








h 2 


cosh(7) 


v / 3cosh(7+^l) 


[— 3ei+2 cosh(7)] cosh(7) 





sinh(7) 


4rsinh(7— 70) 


— ei sinh(7)+sinh(27) 



Table 4: The Hamiltonians couplings for the four branches where 7 = — ei), 70 = 
w = exp (^y 1 ) and ei = e 2 = ±1- The variable r is defined as r = sinh(7) sinh(7 — 70) . 
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The eigenspectrum of the transfer matrix T(X) can be determined by using the algebraic Bethe 
ansatz framework developed in [5] . Here we shall present only the main results that are necessary 
to unveil the thermodynamic limit properties of the spin-1 chain associated to the vertex model 
2B. Its transfer matrix commutes with total spin operator Yli=i &i where SI denotes the azimuthal 
component of a spin-1 matrix acting on the i— th site. Consequently, the Hilbert space can be 
separated in 2L + 1 sectors labeled by an integer number n = L — r where r = 0, ±1, . . . , ±L. We 
shall denote by A„(A) the eigenvalue of the transfer matrix T(A) on a given sector n. By applying 
the method [5] we find that the expression for A n (A) is, 

sinh[Aj — A + 



An(A) = }\< 



12 J 



+ 



+ 



s sinh[A,- - A - ] 

sinh(A) 1 L sinh[2(A - A,) + ifi] sinh[A - A,- - ^ 



sinh (A + *f>) 



2 £l sinh[2(A - A,) - Hf] sinh[A - A, + ^] 



sinh(A + j^) sinh(A) 



-, L 



n ei smh[A - A -- f 5 + ^ ] (150 

f = \ sinh[A - X j + ^f] 



cosh(A - 2^) cosh(A) 
provided that the rapidities {A} satisfy the following Bethe ansatz equations 



We now have the basic ingredients to investigate the ground state and the nature of the 
excitations of the corresponding spin chain. The corresponding eigenvalue E n (L) is obtained by 
taking the logarithmic derivative of A n (A) at A = and using as normalization J = ~ The 
result is, 



n 2 " f 

E " [L) = 61 |r cos(f ) - n col(2A,) < 1K > 



From now on we shall concentrate our attention to study the properties of the anti-ferromagnetic 
regime e± = +1. Further progress can be made after the identification of the distribution of roots 
{Xj} that are able to reproduce the low-lying energies of the spin Hamiltonian. This step is done 
by numerically solving the Bethe ansatz equation (I15ip and substituting the roots in Eq. (I152[) . 
We then compare the results for E n (L) with the exact diagonalization of the Hamiltonian up to 
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L = 12. We observe that though the Hamiltonian is not hermitian we found that its eigenvalues 
are all real. By performing this analysis we find that the ground state for L even sits on the sector 
with zero magnetization n = L. Interesting enough, the shape of the roots on the complex plane 
depends whether |f is even or odd. In Figure [5] we exhibit the Bethe roots for L = 8, 10, 12 and 14 
scaled by the factor -. 
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Figure 5: The ground state roots Zj = ^ for (a) L = 8, (b) L = 10, (c) L = 12, (d) L = 14. 

We then solve the Bethe equations (I15ip for larger values of L to figure out the pattern of the 
roots in the thermodynamic limit L — > oo. In Figure [6] we show the corresponding Bethe roots for 
L = 40 and 42. For better display of the roots curvature we shown the positive and the negative 
imaginary parts separately. 

The above analysis leads us to conclude that as L — > oo the roots Xj cluster in a complex 
2— strings. Each 2— string has the same real part fj, s and equally spaced imaginary parts, 

ITT 

Xj = fij ± — , /Xj e E. (153) 
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Figure 6: The ground state roots Zj = — for L = 40 and L = 42 where (a) Zj = Xj + ii/j and (b) 
z j = x j ~ l Dr 



We now substitute the structure (11531) in the Bethe equations ( 1151ft and by taking its logarithm 
we find that the resulting relations for 



5n 
12 



TV 



n 

-2irQj + ( 2 Vj - 2 Hk, |) , for j = 1, 



n 



(154) 



where the function <j)(x,y) = 2 arctan[tanh(x) cot(y)] and Qj are integer or semi-integer numbers 
characterizing the logarithm branches. 

Considering our numerical analysis we find that the low-lying spectrum can be described by 
the following sequence of Qj numbers, 



Qj 



L 



— n — 1 



j - 1 , for j = 1, 



L 



— n. 



(155) 



For large L, the roots approach toward a continuous distribution with density p(/x) given by, 



(156) 
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where the counting function is 

1 



2n 



(157) 



Strictly in the L — >• oo, Eqs. (I156]I157|) go into a linear integral equation for the density p(fj>) which 
can be solved by standard Fourier transform method. The final result for p(fx) is, 

p(a) = l — — -. (158) 

HK ^ J vrcosh(2 / u) v ; 

The ground state energy per site e ro can now be computed by taking the infinity volume limit of 
Eq. fll52p and with the help of Eq. (11581) . By writing the result in terms of its Fourier transform 
we find, 

sinhjg] - sinhpf] _ _2 V3 
n au cos h[^] sinh[^] ~ vr + 9 ' 1 ' 

Let us turn our attention to the behaviour of the low-lying excitations. These states are 

obtained from the Bethe equations (11541) by introducing vacancies on the sequence of numbers 

Qj. This method is by now standard to integrable models and for technicalities see for instance 

[23] . It turns out that the energy e(/x) and the momenta p(/i), measured from the ground state, 

of a hole excitation are given by, 

POO 

e(ju) = 27rp(/x) and p(/i) = / e(x)dx. (160) 

J ii 

In order to calculate the low-lying dispersion relation e(p) we have first to compute the integral in 
Eq. (ll60p and afterwards eliminate the variable \i from function e(/i). By performing these steps 
we find that the dispersion relation is 

e(p) = 2sin(p). (161) 

From Eq.f ll6ip we conclude that the low momenta excitations have a linear behaviour on the 
respective momenta. Therefore, the excitations are massless corroborating the connection between 
the geometric form of the curve and the nature of the spin-1 chain excitations made at end of 
section HI 
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6 Conclusion 



In conclusion, we have investigated the solutions of the Yang-Baxter equation for nineteen vertex 
models invariant by the action of parity-time reversal symmetry. We have developed a method to 
solve the corresponding functional equations from an algebraic point of view without the need of a 
priori spectral parameterization assumption. The structure of the algebraic manifolds constraining 
the Boltzmann weights follows a rather universal pattern allowing us to classify the possible 
integrable vertex models in four different families. They share the same fundamental algebraic 
curve which turns out to be a conic depending on three basic weights. By using the standard 
parameterization of conies we are able to obtain the dependence of the Boltzmann weights on the 
spectral parameter from a unified perspective. Three of such vertex models were already known 
before but the fourth model appears to be new in the literature. 

We have observed an intriguing relation between the form of the main algebraic curve and 
the nature of the low-lying excitations of the related spin-1 chains. In the regime in which the 
geometric form of such curve is an ellipse the excitations are massless while when we have a 
hyperbolic structure the excitations have a mass gap. This fact is supported by previous knowledge 
on the physical properties of the first three spin-1 chains and also by the exact solution of the 
novel nineteen vertex model. 

It seems interesting to investigate whether or not the scenario described above remains valid 
for other families of VT invariant vertex models with larger number of states. In principle, 
the systematic method developed here can be extended to tackle other models whose statistical 
configurations preserves at least one £7(1) symmetry. Another interesting problem is to study the 
algebraic invariants underlying vertex models that are not invariant by the VT symmetry and in 
particular to unveil the form of their principal algebraic curve. We hope to address these questions 
in future publications. 
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Appendix A: Three terms relations 



This appendix is devoted to the presentation of some technical details entering the solution of the 
group G^± relations. An immediate consequence of the last invariant (1511) is that the weight 
can be easily written either in terms of the weights a+ and b^} or by means of the weights a$ 
and b_ , namely 



r ± &<*> 



(A ± - r±) ag - 



for % = 1,2. 



(A.l) 



It turns out that Eq.f lA.ll) can be used to obtain a similar expression for the weig ht d (i \ In 
fact, by substituting Eq.( 1A.lj) in Eqs.( l26K9ll50l) and by carrying on some simplification one finds, 



(i-A±r± + r|) 





to 


• 


(*)" 


'bf 


1 




<4 



-A ± a®b® + 



/Al 2 

■ w 



;-A ± + r ± )a« + ^ )_ 



for i = 1,2. 



(A.2) 



By using the form of the invariants A± given by Eq. fl22j) in the expression (1A.2j) one is able to 

i 2 



take the square root of the weight [dW] . The final result is 



i - A ± r ± + r| 



A H 



-A ± + r±) + 6± 



for i = 1,2. 



(A.3) 



This means that the compatibility among the invariants T± is equivalent to the matching of the 
expressions for the weights /w and dy> coming from the distinct charge sectors ±. By imposing 
such consistency one finds that Eq. llA.lj) requires us to solve, 



(A_ - r_) a W - 6 W 



a? - I\ 6 



(A + -T + )a + 



« 6 w 



,(*) 



r_6 



o 



for % = 1,2 (A.4) 



while Eq.( 1A.3l) implies that the compatibility for rfW is, 



A_(i-A + r + + q 



a + (i - A_r_ + r 2 _) 



-A_ + T_) a<? + &<? 
-A + + T + )a i l ) + b { l ] 



fe » c « _ q for j = lj2 . (A.5) 
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The analysis of Eqs. (1A.4)A.5P for the index i = 1 is immediate since the weights a_ and c_ are 
still free. Therefore, Eqs.( lA.4llA.5l) for i — 1 are easily solved by eliminating the weights and 
c_ , namely 



r_(A + - T + )aJ ) 6 (1) - a^b^ - T.b^b^ + T 



(%(!) 



(A+ - r+)aV 1) b (1) - (A_ - r_)aV iJ ^ " b -' b + + r +( A - " r 



(A.6) 



and 



• VA-(i - A_r_ + r 2 )(i - A+r+ + r 2 + )b { y b 



2u(lk(l) c W 



<{ (A+ - T + )a^b W - (A_ - r_)aV iJ 6V iJ - 6^ + T + (A_ - T 



(A.7) 



On the other hand the solution of Eqs. (lA.4fA.5p for the index % = 2 is more involving. This 

(2) (2) (2) 

is the case because the weights a_ , 6_ and c_ have already been determined in terms of the 
amplitudes a+ , b+\ c+ \ and S 2 ^ thanks to the previous relations (I32ll46p . In addiction to 



that we also recall that the weights and are also determined in terms of a ? , b^ and c 



.(2) 



through Eqs. (1A.lllA.3p . Considering all these information together with the fact that the weight 



.(2) 



-i 2 



can be eliminated with the help of the invariant A + (|22l) one concludes that Eqs. (1A.4IIA.5P 



are in fact a polynomial relation on the remaining weights ao/ and b y + . It turns out the expression 



',(2) 



of this polynomial associated to Eq. (1A.4[) . after using the constraints (I35|45p . is 



Bi 



,(2)' 



■,(2) 



,(2)' 



6? + B A 



■(2) 



(A.8) 



where its coefficients are given only in terms of the invariants values by, 

Bi = ^[A 2 ,r_(-A + + r + ) + A + ^(A + r + + r_r+-r 2 ^ + A + fi + r_fi-r + fi) 
- ^(r^ + A + fi + 2r + fi + 2fi 2 )], 



(A.9) 



b 2 = A^r+(A + + 2r_ -r+) + A^(i-2A+r_ - A + r + + r_r+-4r + fi) + A+^ 2 (-i 

- A+r_ + A+r+ - r_r + + A^r_r + + A+n + r_fi - r + n - A 2 + r + n - A + r_r + ^ 
+ A + r 2 + n - o 2 - A+r+fi 2 - 2r_r + fi 2 + r 2 + n 2 ) - ^ 3 (i - 2A + r + + r 2 + + a 2 + t 2 + 

- 2T + tt + A + r 2 fi + 2fi 2 -4A + r + fi 2 -4r + fi 3 ), (A.10) 
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b 3 = -A3(A + + r_-r+) + A^(r_ + A^r_-r+-A + r_r + + 2fi) + A + ^ 2 (A + + 2r_ 
- 2r + - A^r+ - 2A + r_r + + A+r^ - ft - r_r+ft + A+ft 2 + r_ft 2 - r+ft 2 ) 
+ ^^A+r 2 ; - 3ft + 4A+r+ft + r 2 h ft- A+ft 2 + 2r+ft 2 - 2ft 3 ), (A.n) 

b a = -A 3 f r + (i + r_r + ) + A^(r_-r + + A + r 2 ; + 2r 2 h fi) + A + ^ 2 (r_- A + r_r + 
+ t_t\ - ft + A + r + ft - r_r + ft + A + r_r^ft + r + ft 2 + r_r^ft 2 ) - ^ 3 (ft - 2A + r+ft 
+ r 2 K ft + A 2 ^r 2 ;ft^ 3 -2r + ft 2 ^ 3 + 3A + r 2 ;ft 2 ^ 3 + 2r 2 K ft 3 ^ 3 ). (A.12) 

(2) (2) 

As argued in the main text we are searching for solutions in which a + and b + are independent 
of each other. This means that we have to set all the coefficients of the polynomial flA.81) to zero. 
By imposing this condition for B\ one is able to eliminate the invariant T_ whose expression can 
be further simplified with the help of the constraint ( |36|) . The final result is, 

r- = ^-(r + + ft). (A.13) 



By substituting the result ( 1A.13j) in the expressions of the remaining coefficients (lA.10HA.T2l) 



and by once again taking into account the constraint (I36I) we find that B2, -B3 and B4 vanish. The 
same reasoning can be repeated for Eq. (lA.5p when i = 2. We conclude that it does not impose 
additional restriction besides Eq. (IA.13~j) . 

We conclude by observing that the weights a®, b_ and can also be written in terms of 
the amplitudes a+, b+ and c+. In order to see that we first substitute the weights /W e 
considering the index + of Eqs. (IA.4|A.5p . in the expressions of the weights b_\ b_ , c_ , " l2) u 



a_ , a: 



and c^, see Eqs. (l30|l32ll46IIA.6IIA.7l) . We next use the constraint A + (|22|) to eliminate the amplitude 



^ 2 



besides the explicit form of the relations between the invariants (I27f35ll36f45f A. 131) . By 
performing such steps we are able to write rather simple expressions for the weights o_ , b^\ 
for i = 1,2. They have given in the main text, see Eqs. fl56H58l) 



39 



Appendix B: Particular invariant solutions 



We shall describe particular solutions for the invariants coming from branch 1A. Besides \I/ = 
A + = 1 we have two other possibilities that provide us non-null Boltzmann weights. We shall 
denominate such special branches as follows, 

• Branch IS : * = -A+ = 1, (B.l) 

• Branch 2S : * = A+ = -1. (B.2) 

It turns out that the choices (1B.l|B.2j) lead us to the situation in which we do not have any 
free parameter at our disposal. In Table |6] we present the corresponding invariant values. 

By comparing Table |6] with Table |2] it is not difficult to see that the branch IS is a particular 
case of branch IB once we choose A + = ±v3. By the same token comparison with Table [3] reveals 
us that branch 2S becomes equivalent to branch 2 A provided we set A + = ±2. 
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